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Abstract —The economic dispatch of wind power units is quite 
different from that in conventional thermal units; since the 
adopted model should take into consideration the intermittency 
nature of wind speed as well. Therefore, this paper uses a 
model that takes into account the aforementioned consideration 
in addition to whether the utility owns wind turbines or not. 
The economic dispatch is solved hy using one of the modern 
optimization algorithms: the particle swarm optimization 
algorithm. A 6-hus system is used and it includes wind-powered 
generators besides to thermal generators. The thorough analysis 
of the results is also provided. 
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I. Introduction 

As a result of the intermittency nature of wind power, 
the economic dispatch of resources that include wind power 
is quite different of that in pure conventional thermal units. 
Therefore, how can the economic dispatch of this promised 
future resource of energy be achieved? 

Various mathematical programming approaches used to solve 
this kind of optimization problem in power systems, based on 
linear and nonlinear programming were proposed, including 
Newton method, quadratic programming, and interior-point 
method 111. The mathematical methods utilize the first or 
second derivative information in essence. In this way, it is 
apt to fall into local optima. Eurthermore, there is a difficulty 
of applying gradient-based optimization techniques. Therefore, 
various non-classical optimization methods have emerged to 
cope with some of the traditional optimization algorithms’ 
shortcomings. The main modern optimization techniques are 
genetic algorithm (GA), evolutionary programming (EP), ar¬ 
tificial neural network (ANN), simulated annealing (SA), ant 
colony optimization (AGO), and particle swarm optimization 
(PSO). They have been successfully applied to wide range 
of optimization problems in which global solutions are more 
preferred than local ones |2l,|3l. 

Kennedy and Eberhart first introduced particle swarm op¬ 
timization (PSO) in 1995 as a new heuristic method 0|. 
In |3l there is a comprehensive coverage of different PSO 
applications in solving optimization problems in the area of 
electric power systems up to 2006. The review in Q is about 
the historical research production of the economic dispatch 
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considering the wind power, besides that it also presents some 
models and different optimization algorithms as well. 

In 2008 m is one of the pioneer studies about the economic 
dispatch including the wind power was reported. It also 
includes the definitions about the wind power cost and its 
factors in wind energy conversion systems (WECS) combining 
both cases, whether the operator owns WECS or not. In 
addition to the direct cost of wind power, cost factors of the 
overestimation and underestimation of wind power have also 
been proposed. 

This paper is intending to investigate the interconnection 
of wind generators besides the conventional generators into 
power systems and its impact on the generation resource 
management. The rest of this paper is organized as follows. 
Section El introduces the motivation and the problem state¬ 
ment. Section mil it discusses the analysis and characterization 
of wind speed and power. Section |IV] the particle swarm 
optimization (PSO) algorithm is described. In Section |V] the 
results of the implementation of PSO to find the economic 
dispatch of a benchmark system are discussed to some extent. 
Einally, Section IVTl gives the conclusions. 

II. Problem Eormulation 
A. Problem objectives 

The objective function of the optimization problem in this 
paper is to minimize the operating cost of power generation 
power from a combination of wind-powered and conventional 
generators. The operating cost of conventional thermal gener¬ 
ators is represented by a quadratic equation as following Q: 

Ci — biPi Ci ( 1 ) 

Where p; is the generation power from the ith conventional 
generator; and a, b and c are the operating cost coefficients of 
the Ith generator. The wind power generation cost which 
may be not exist if the power operator owns the wind powered- 
generators, but it could be considered as a payback cost or a 
maintenance cost 0: 

— diWi ( 2 ) 

Since w; is the scheduled wind power from the ith wind- 
powered generator; and d; is the direct cost coefficient for 
the Ith wind generator. 


Because of the uncertainty of generated wind power, there 
are two scenarios of wind power costs. The surplus of wind 
power as a result of underestimation of the available wind 
power and hence scheduling the wind power Wi less than what 
it would be. Thus, Cp appears as a penalty cost Ih). 

pWr,i 

Cp^i = / {w- Wi)fyj{w)dw ( 3 ) 

J Wi 

Where is the Weibull distribution function for wind power, 
for more details see section (Illl-Cb : and kp is the penalty cost 
coefficient. 

On the other hand, a deficit of wind power which occurs by 
the overestimation of the available wind power and scheduling 
the wind power Wi more than it would be available. At that 
situation, the deficit will be compensated by a reserve power 
sources. That means there is also a cost for the deficit of wind 
power. Thus, Cr is presented as a reserve cost. 

pWi 

Cr,i = K^i / {Wi - w)fn,(w)dw ( 4 ) 

^0 

kr^i is the reserve cost coefficient for the ith wind generator. 
B. Problem Constraints 

Due to the physical or operational limits in practical sys¬ 
tems, there is a set of constraints that should be satisfied 
throughout the system operations for a feasible solution [jS). 

• Generation capacity constraints: 

For normal system operations, real power output of each 
generator is restricted by lower and upper limits as 
follows: 

pT"" < k < pT°''' 

0 < Wi < Wr,i 

Since Wr^i is the rating wind power from the ith wind- 
powered generator. 

• The transmission line losses constraints: 

Q ^ Qmax 

^line,i ^ ^Une,i ID 

Siine,i is losses of the Ith transmission line. 

• Power balance constraint: 

The total power from conventional and wind generators 
must cover the total demand. 

M N 

'^Pz + '^Wi = D ( 8 ) 

i=l i=l 

Where M number of conventional power generators; N 
number of wind-powered generators; and D is the demand 
which equals to the system load and losses. 


(5) 

(6) 


C. Problem Statement 

In summary, the objectives of optimal economical dispatch 
is to minimize the operating cost from the conventional and 
wind-powered generators includes the penalty of underestima¬ 
tion and overestimation of wind power, subject to the certain 
constraints. 

The model of economic dispatch for thermal and wind- 
powered generators ||6]: 

M N N N 

Ci{pi) + '^Cw,i{Wi) + '^Cp^i{Wi) + '^Cr,i{Wi) (9) 

i — 1 2=1 2 = 1 2 = 1 

subject to: The constraints that are represented as in equations 

© - ®. 

Note using a classic economic dispatch approach for the 
model in equation ©, which takes the partial derivative of the 
objective function respect to generator outputs; it’s difficult 
due to the integrals in the wind power cost terms as in 
equations © and @, Therefore, Particle Swarm Optimization 
(PSO) algorithm is used for solving this optimization problem. 

III. The Analysis of Wind Speed and Power 
A. Probability Analysis of Wind Power 

Before starting the discussion of economic dispatch of 
systems that contain wind-powered generators, it will be a 
good idea to identify the wind speed characterization by 
probability principles and its subsequent transformation to 
wind power. 


B. Wind Speed Characterization 

The wind speeds in a particular place can be considered as 
a Weibull distribution over time The probability density 
function (pdf) of the Weibull distribution fv{v) is given by: 



0 < u < c» (10) 


Where v is the wind speed; c is scale factor; k is the shape 
factor. Fig. [T] illustrates the Weibull pdf with shape factors 
k=2, and curves of scale factor c = 5 m/s, 15 m/s, and 25 
m/s are indicated. 



Fig. 1. Weibull pdf of wind speed for several values of scale factor c (9) 








The cumulative distribution function (cdf) of Weibulll distri¬ 
bution Fv{v) is obtained by integration of pdf: 

pV 

Fv{v) = / fv{v) dv = I — (11) 

Jo 

C. WECS Input/Output and Probability Functions 

For wind energy conversion systems (WECS) as it is shown 
in Fig. |2] the wind power curve from probability point of view 
can be represented in three regions as in equation (fTSl i 161 . 

{ 0; {v < Vi or V > Vo) 

Wr ; iVi<V< Vr) (12) 

Wr\ [Vr <V < Vo) 

Where w is the wind power; Wr is the rating power of WECS; 
Vi is the cut-in wind speed; Vo is the cut-out of wind speed; 
Vr is the rating wind speed at which the rating power Wr is 
captured. 



Fig. 2. The captured wind power curve 


The linear transformation from wind speed to wind power 
in the linear region {vi < v < Vr) is done as following ||6l: 

rri \ t ^ {w-b) 

. w = 1 (V) = av + b . .V = 1 (w) => u =- 

a 

So now V in terms of w, and thus fv{v) =i> fw{w) as follows: 


.-. fw{w) = fv 


( w — b\ 

1 

V a j 

a 


(13) 


where: 

T is the general transformation; w wind power random vari¬ 
able; V wind speed random variable; 

For Weibull distribution function, the transformation will 
lead to discrete and continuous ranges as following: For 
discrete portions: 


Pr{W = 0} = Fviv,) + {l-Fvivo)) = 

^(14) 

Pr{W = Wr} = Fv{vo) — Fv{vr) = 

(15) 

While for the continuous portion of the wind power curve: 


fw{u>) 


klv^ f {1 + pl)Vi V^~^'^ 

WrC \ C J 


(16) 


w 



^ _ (Vr - Vi) 
Vi 


IV. Particle Swarm Optimization Algorithm 

The PSO is originally suggested by Kennedy and Eberhart 
based on the analogy of swarm of bird and school of fish 
0. The algorithm was simplified and used for solving the 
optimization problems. 

A. Standard PSO Algorithm 

The following is the conventional terminology of the pa¬ 
rameters in PSO: Let x and v denote a particle coordinates 
(position) and its corresponding speed magnitude (velocity) 
in a search space, respectively. Therefore, the ith particle is 
represented as Xi = [xii,Xi 2 , ■■■■,Xim]- Since m is the last 
dimension or coordinate of the position of the the ith particle 
in the search space and so that the dimension = 1, 2,..., m. 
The best previous position of the ith particle is saved and 
represented as m, 

pbesU = [pbesUi,pbesti 2 , ..■.,pbesUm]- 
The position of the best particle among all the particles in the 
group is represented by the gbest. In a particular dimension d 
there is a group best position which is gbestd- 
The velocity for the ith particle is represented as, Vi = 
\vii,Vi 2 , ■■■■,Vid]- The modified velocity and the position of 
each particle can be calculated by using the following formu¬ 
las: 

= w*v^d + ci*U*{pbest'lj^-x'lji) + C2*U*{gbest'l-x^,j) 

(17) 

(18) 

z = 1, 2,...., n; d=l,2,...,m 

Where x^^,vf^ the position and the velocity of the ith particle 
in the </th dimension at an iteration k; n number of particles in 
a group; m number of members in a particle; w inertia weight 
factor; ci,C 2 acceleration factors; U uniform random number 
in the range [0,1]; 

The velocity should between if y^ax 

is too high, particles might move past good solutions. While 
if 11 ™“^ is too small, particles may not explore sufficiently 
beyond local solutions. In many experiences with PSO, was 
often set at 10 - 20% of the dynamic range of the variable on 
each dimension Eol. 

The constants ci and C 2 represent the weighting of the 
stochastic acceleration terms that pull each particle toward the 
pbest and gbest positions. The acceleration constants ci and 
C 2 are often set to be 2 according to past experiences ifTol . 
Suitable selection of inertia weight w in equation (ITtT i provides 
a balance between global and local explorations, to find a 
sufficiently optimal solution. In general, the inertia weight is 
set according to the following equation: 

(v)inax VJfoin) 

W = Wmax -T- X iteT (19) 

'iteVmax 

Where iteXmax is the maximum number of iterations (gener¬ 
ations). PSO algorithm flowchart is shown in Fig. [3 

















Fig. 3. PSO algorithm flowchart 


V. PSO Implementation on IEEE Benchmark System 

PSO algorithm is coded in MATLAB to determine the 
economic dispatch of the generation power for an electrical 
power system with considering various levels of wind power 
penetration. In this paper, the PSO algorithm is implemented 
for solving the economic dispatch of a 6-bus system that 
includes wind-powered generators in addition to the con¬ 
ventional generators. Since the used model in equation (|9]l 
has an objective function that contains some terms that need 
integration, the PSO algorithm can be applied efficiently to 
solve this type of objective functions. 

A. The Data of The System 

The 6-Bus System in Pig. |4] is adopted to calculate the 
optimal economic dispatch of generation power including wind 
power. This system consists of six buses and four generators 
at buses 1, 2, 3, and 4, generators at buses 3 and 4 are wind- 
powered generators. There are seven transmission lines, and 
there are no ETC transformers or VAR compensation devices 
in this system. Por all detailed data of this 6-bus system refer 
to HD. 



Fig. 4. Single-line diagram of 6-bus system (m 


As TABLE U is shown, the system has two conventional 
thermal units and two wind-powered generators. The higher 


TABLE I 

Generators data of 6-bus system 


Gen. No. 

a 

($/MW'‘2.hr) 

b 

(S/MW.hr) 

c 

Pg low 
(MW> 

Pg high 
(MW> 

1 

0.012 

12 

105 

50 

250 

2 

0.0096 

9.6 

96 

50 

250 

3 

0 

8 

0 

0 

40 

4 

0 

6 

0 

0 

40 


output of each wind-powered generator is 40 MW. While 
the direct cost of wind power 8 and 6 %/{MW.hr) for 
wind-powered generator 3 and 4 respectively. The difference 
in direct cost of wind power is for including the variety 
purpose in the study, and to get more options of dispatching. 
The parameters of the wind turbine are cut-in wind speed 
Vi = 5m/s, rating wind speed Vr = 15m/s, and cut-out wind 
speed Vo = 4:5m/s. 

B. The Objective Function 

The aim is to implement the PSO algorithm with the model 
which is discussed in section HI] The objective function is 
presented by equation ®, the term '^^e cost 

of the real power of thermal-generators, J2iLiC^w,i{wi) is 
the direct cost of wind power, '^^e penalty 

cost of the underestimation of the available wind power, and 
Sfci Cr,i{wi) is the reserve cost of the overestimation of the 
available wind power. The two latter terms in the objective 
function have integrals as they are represented in equations 
01 and 0. 

Pig. |5] illustrates the cumulative probability distribution of 
wind power, it is produced from the integration of equation 
[Thl This figure has such importance for the next investigations. 
Note; w is the available wind power, while Wr is the rating 
wind power. 



Fig. 5. Cumulative probability distribution of wind power vs. normalized 
wind power 

C. PSO Solution for Base Case 

The base case with total load of 400 MW. The parameters 
of Weibull distribution of wind speed here are scale factor 
c =5m/s, while the shape factor k is 2. By assumption, 
the reserve cost coefficient as a result of overestimation of 
available wind power would be 1 %/MW.hr. On the other 
hand, penalty cost coefficient as a result of underestimation of 
available wind power is 0 %/MW.hr, this means the utility 
owns wind turbines so there is no penalty of surplus produced 
wind power. The changing of these coefficients and their effect 
on the total cost will be investigated later. 





























The economic dispatch of the base case is in TABLE [H] 
The minimum cost of real power from both thermal and wind 
power generators is 4777.49 %/hr. Wind-powered generators 
in these conditions supply maximum outputs because they 
are more economic. While the first thermal generator supplies 
less power to the system than the second generator because 
its generated power is more expensive. All generators’ cost 
coffeceitns are in TABLE |I] 

As it is shown in TABLE|II]the outputs of generators are equal 
to the demand plus losses in the system. 

TABLE II 

PSO RESULT OF ECONOMIC DISPATCH FOR BASE CASE (TOOMIP) 


( 3 ) kf — 1 


300r 



(b) kr=lO 


25 


O 150 



•-G1 
—G2 
—G3 
-G4 


C Weibull scale factor 


25 


Pgi 

Pg2 

Pg3 

PG4 

Losses 

Cost 

(MW) 

(MW) 

(MW) 

(MW) 

(MW) 

($/hr) 

92.82 

230.64 

40 

40 

3.462 

4777.49 


C =5 m/s. K=2. (Weibull PDF parameters) 

Ki=l S/MW.hr. Kp=0 S/MW.hr as utility owns wind turbines (wind power cost factors). 


D. PSO Solution for Different Loading 

The economic dispatch solution by PSO algorithm when the 
system load increases gradually is as in TABLE |III] The limits 
of the transmission lines losses are neglected because they are 
relatively low. 

TABLE III 

PSO RESULT FOR ED OF DIFFERENT LOAD CASES OF 6-BUS SYSTEM 


Load 

Pgi 

Pg2 

Pg3 

Pg4 

Cost 

(MW) 

(MW) 

(MW) 

(MW) 

(MW) 

(SAD 

400 

86.66 

233.34 

40 

40 

4728.225 

450 

120 

250 

40 

40 

5448.225 

500 

170 

250 

40 

40 

6222.225 

C=5ni/s, K=2. (Weibull PDF parameters) 

K,=l S/MWJir, K«=0 $/MWJii. 


E. The Effects of Wind Power Cost Coefficients 

Next the variety of wind power cost coefficients and wind 
speed factors will be investigated. So that their effects on the 
output schedule of the generators and hence the total cost for 
the base case 400 MW are presented as following. 

The shape factor of wind speed probability Weibull distribu¬ 
tion k = 2 and it is kept constant at this value. While the scale 
factor c is changing between 5m/s to 25 m/s. The constant 
direct costs of wind power from wind-powered generators 
3 and 4 are 8 and 6 %/(MW.hr) respectively. Eor sake of 
convenience, hereinafter the units will be dropped from these 
coefficients. 

1) The Effects of Reserve Cost Coefficient: Eirst, assume 
that the utility owns wind turbines, so that the penalty cost of 
additional available wind power over scheduled power will be 
0, and this also means the coefficient fcj, = 0 as it is derived 
from equation Q. 

Fig. |6] shows the result of PSO algorithm for the economic 
dispatch of generators’ outputs as a function of the scale factor 
of Weibull distribution of wind speed c for different values of 
reserve cost coefficient kr. 

When c scale factor of Weibull distribution of wind speed 
increases, the reserve cost decreases. That can be verified 
from Pig. |5] When c increases, the probability of wind power 
decreases, then the reserve cost reduces as well. Thus, the 


(c) fcr = 100 

Fig. 6. Generators’ outputs vs. Weibull scale factor (c) for some values of 
reserve cost in base case 400 MW 

reserve cost reduces by increasing c, the scheduled outputs 
of wind-powered generators will increase gradually as in 
Pig. |6bl In Pig. |6c] there is a small increase in wind-powered 
generators form c =20 to c =25, it is a small change because 
the reserve cost coefficient in this case relatively high kr=lQQ. 

2) Critical Reserve Cost Coefficient: Pig. |7] shows the 
outputs of generators for a variation of reserve cost coefficient 
kr for two values of scale factor c = 5 and c = 20 in order 
to see where the critical change in wind power schedule 
begins. In Pig. [Ta] when c=5, the critical change in wind 


1^250 
1^200 
o )50 
IfOO 


5 10 15 20 25 30 35 40 45 50 

Kr Reserve cost factor, C=20 

(a) c=5 (b) c=20 

Fig. 7. Generator outputs vs. reserve cost coefficient (kr) for two values 
of scale factor (c) in base case 400 MW 

power schedule starts when kr = 6 for generator (3) and 
kr = 8.6 for generator (4). The drop of outputs happens in 
generators (3) before generator (4) because generator (3) has 
a higher direct cost (8 %/MW.hr) than that of generator (4) 
(6 %/MW.hr). While in the other case when c=20 as in Pig. 
|7bl the change of wind power scheduling occurs at higher 
values of kr because in this case, the scale factor c of Weibull 
distribution of wind speed has a higher value. 



3) The Effects of Penalty Cost Coefficient: When kr = 0 
and kp f- 0, the schedule of generators as in Pig. [8a] for 
various values of kp remains constant for different values of 
scale factor c. In this case it should get all available wind 
power since there is a penalty cost for a surplus wind power. 
Pig. [8b| shows that penalty cost Cp=0, thereby all available 
wind power is scheduled from both wind-powered generators 
as Pig, nil shows. 


4) The Effects of The Reserve and Penalty Cost Coefficients: 
The effect of both of the reserve and the penalty cost coef¬ 
ficients of the economic dispatch outcomes) is illustrated in 
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Kp Penalty cost factor, When C=5 or C=20 
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Kp Penalty cost factor, when C=5 or C=20 


(a) Generators’ Outputs (b) Penalty Cost Cp 

Fig. 8. Generators’ outputs and penalty cost Cp vs. penalty cost coefficient 
kp for two values of scale factor c in base case AOOMW 


Fig. |9] for the base case with a scale factor c=5 and the wind 
cost coefficients are not equal to zero (kr ^ 0 and 7 ^ 0 ). 
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Kp Penalty cost factor, C=5 
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Kp Penalty cost factor, C=5 


electric power operators. The used model of the economic 
dispatch that consideres the wind power uses the manipulation 
of the probability of the underestimation and overestimation 
of the availability of the wind power. It also takes into account 
whether the utility owns wind turbines or not; these are the 
main features of this model. The variations of wind speed 
parameters and their impacts on the total cost investigated by 
6 -bus system, some valuable conclusions have been noticed. 
The incremental reserve and penalty costs of the available 
wind power can be compared to the incremental cost in the 
conventional-thermal units that have a quadratic cost; this com¬ 
parison could lead to useful simplihcations of the economic 
dispatch models that include thermal and wind power. PSO 
algorithm needs some work on selecting proper parameters 
and it also needs some further mathematical description for 
its convergence. 

References 


(a) kr = 0 (b) kr = 60 

Eig. 9. Generators’ outputs vs. penalty cost coefficient kp for some values 
of reserve cost coefficient kr in the base case AOOMW and c = 5m/s 

When the reserve cost coefficient kr increases, the sched¬ 
uled wind power decreases. Until no scheduled power comes 
from wind-powered units when kr > 60; because the high 
value of kr makes the wind power to be not an economic 
option. Thereby, all the scheduled power comes from thermal- 
powered units for any value of kp, as in Fig. |9b] 

Fig. [TOl can be considered as a part of Fig. |9] when kr=2Q 
but now for two higher values of the scale factor c =10 and 
c= 20 . 
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(a) c = 10 (h) c= 20 

Fig. 10. Generators’ outputs vs. penalty cost coefficient kp for two values 
of scale factor c when kr =20 


As it is shown in Fig. [TOl with a higher scale factor c of the 
probability distribution of wind speed, the outputs of wind- 
powered generators become higher as well. Furthermore, Fig. 
llOal illustrates that the wind power outputs will increase with 
higher values of the penalty cost coefficient kp. 

This is what happens when the utility does not own the wind 
turbine, therefore the scheduled wind power is produced as a 
compromise between the penalty cost and the reserve cost of 
wind power. 


VI. Conclusion 

The PSO algorithm is appropriate for hnding the optimal 
economic dispatch of generators’ outputs that include wind- 
powered generators. Moreover, PSO algorithm is simple and 
easy to apply since it has fewer number of parameters to 
deal with comparing to other modern optimization algorithms. 
PSO can be applied in wind power bid marketing between 
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